Statistical theory of a quantum emitter strongly coupled to Anderson-localized modes 
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A statistical theory of the coupling between a quantum emitter and Anderson-localized cavity 
modes is presented based on a dyadic Green's function formalism. The probability of achieving the 
strong light-matter coupling regime is extracted for an experimentally realistic system composed 
of In As quantum dots embedded in a disordered photonic crystal waveguide. We demonstrate 
that by engineering the relevant parameters that define the quality of light confinement, i.e. the 
light localization length and the loss length, strong coupling between a single quantum dot and an 
Anderson- localized cavity is within experimental reach. As a consequence of disorder-induced light 
confinement provides a novel platform for quantum electrodynamics experiments. 
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The ability of controlling the coupling between a sin- 
gle emitter and a photon grants access to a rich set of 
quantum optical phenomena [l| and provides a route to- 
wards quantum networks for quantum information tech- 
nology 0. The coupling strength can be enhanced by 
placing the emitter in a small optical cavity that traps 
the emitted photon jsj. As a result, the spontaneous 
emission rate is enhanced through the Purcell effect and 
if the coupling strength exceeds the cavity loss rate, the 
strong coupling regime is reached where coherent oscilla- 
tions between photon and emitter induce quantum entan- 
glement [4|. In solid state systems, strong coupling has 
been achieved, e.g., by embedding semiconductor quan- 
tum dots in highly engineered photonic crystal cavities 
with small mode volumes and high Q-factors |5|. How- 
ever, the fabrication of such cavities requires nanoscale 
accuracy and is therefore inherently sensitive to unavoid- 
able disorder introduced in the fabrication process. 

An alternative route to efficient light confinement 
based on disordered photonic nanostructures has re- 
cently proven to be very promising. In such systems 
the light propagation is strongly modified by multi- 
ple scattering, which can lead to the spontaneous for- 
mation of Anderson-localized modes. These modes 
are robust against fabrication imperfections and prevail 
even after averaging over all configurations of disorder 
[^.Disordered photonic crystal waveguides are well suited 
for investigating Anderson localization since the light 
propagation is primarily one-dimensional (ID) and lo- 
calized modes are formed when the sample length is suf- 
ficiently long. Impressive Q-factors as high as 600.000 
have been reported in such disorder-induced modes 0, 
that subsequently were proven to be due to Anderson 
localization by carrying out the proper configuration av- 
erage 

Recently, Anderson-localized modes were shown to en- 
able very pronounced cavity quantum electrodynamics 
effects by probing the dynamics of single quantum dots 
embedded in disordered photonic crystal waveguides [To| • 
Here we present a theoretical model for photon-matter 



interaction of a single quantum emitter in a ID disor- 
dered medium in the Anderson-localized regime. In a 
multiple scattering medium, light propagation and emis- 
sion is determined by a statistical process giving rise to 
a distribution of photon- matter coupling rates. Two uni- 
versal parameters fully characterize the propagation of 
light in a ID random medium: the localization length 
scaled to the sample length £/ L and the loss length scaled 



to localization length l/£ [y, [llj, while the emission is 
non-universal in the sense that it also depends on the lo- 
cal correlation length of the medium in the presence of 
the emitter [12, Qjl • We derive a distribution of optical 
modes in the Anderson-localized regime of a ID disor- 
dered medium and based on a dyadic Green's function 
formalism extract the associated distribution of photon- 
matter coupling coefficients. From this model we calcu- 
late the probability of strong coupling of a single quan- 
tum dot to an Anderson-localized cavity mode and the 
dependence on localization length and loss length. Con- 
sidering realistic parameters obtained in disordered pho- 
tonic crystal waveguides [14] we conclude that the ex- 
perimental observation of strong coupling of a quantum 
dot to an Anderson-localized mode is within experimen- 
tal reach. 

In a disordered Wl photonic crystal (PhC) waveguide 
membrane light propagation is effectively ID since it is 
confined to the membrane by total internal reflection and 
to the waveguide by the 2D PhC band gap. Disorder in- 
troduces multiple random scattering of the Bloch waves 
and the localization length, £, denotes the average dis- 
tance between scattering events. The criterion for Ander- 
son localization is £ < L leading to an exponential decay 
of the light intensity emitted by a source after ensemble 
averaging over all realizations of disorder [l5[ • Examples 
of calculated electric field profiles of Anderson-localized 
modes in a PhC waveguide are shown in Fig. [TJa) where 
disorder is introduced in the waveguide by random dis- 
placements of the hole positions in the three rows on 
each side of the waveguide. Anderson-localized modes 
are found to appear only in a narrow frequency interval 
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FIG. 1. (color online), (a) 2D finite- difference time-domain 
calculations of the ^/-component of the electromagnetic field 
(color scale) in a disordered photonic crystal waveguide. We 
have used realistic parameters for a GaAs photonic crystal 
waveguide with lattice parameter a = 260 nm, air-hole radius 
r = 0.29a, and a refractive index of n = 2.76, the latter being 
the effective refractive index of a 150 nm thick GaAs mem- 
brane. Disorder is introduced by displacing the position of 
the holes in the three rows on each side of the waveguide with 
a Gaussian random distribution with standard deviation of 
a = 0.03a. Inset: Closeup view of the geometry of the disor- 
dered waveguide (b) Local density of state (LDOS) for a single 
realization of disorder using a stack of layers with refractive 
indices homogenously distributed in the interval 3.45 ±0.70 
leading to £ = 15 |^m for L = 100 \±m, corresponding to real- 
istic values from experiments on PhC waveguides |14 | 



near the cutoff of the waveguide mode where the density 
of states is high 0, [Tg|, [ItJ • The localization length, loss 
length, and sample length are three essential parameters 
of a ID multiple scattering system that enter into the 
model for light-matter coupling. In the following we will 
use physically realistic parameters and estimate the cou- 
pling of single quantum dot emitters to disordered PhC 
waveguides. 

The local light-matter coupling strength is determined 
by the local density of optical states (LDOS) that is cal- 
culated from the imaginary part of the frequency de- 
pendent dyadic Green's function [L8|. The disordered 
waveguide consists of an L = 100 |Ltm long stack of lay- 
ers with randomly varying refractive indices along the 
^-direction with an average refractive index of n = 3.45 
and a layer thickness of 10 nm. We consider the situation 
relevant for cavity QED experiments where the LDOS is 
dominated by localized quasi-modes and therefore can 
be expressed as p(r,u) = \f(x, y,uj)\ 2 p(z,uj)/A e ff(uj), 



where p(z,ui) accounts for the spatial variation of the 
LDOS along the z direction that is obtained from a self- 
consistent transfer-matrix method [l9j], and uj is the op- 
tical frequency. The transverse effective area A e ^(u) = 
J dxdy h 2 (x, y)\f(x, y, uj)\ 2 is calculated from an ordered 
photonic crystal waveguide, where h(x,y) is the refrac- 
tive index and f(x,y, uj) is the simulated transverse elec- 
tric field profile [20| , and both have been scaled to unity. 
The electric field profile is obtained from 3D band struc- 
ture calculations and the weak dependence on z is elim- 
inated by averaging over one unit cell. All parameters 
used in the calculation of the PhC are listed in Fig. [TJ 
and only the size of the refractive index fluctuations in 
the ID multiple scattering model is varied giving rise to 
different localization lengths £ = 7.40/(An) 2 |Ltm for the 
L = 100 |Ltm sample length and the particular form of 
disorder introduced in the present analysis. 

An example of the spatial and spectral variations of the 
LDOS for a single configuration of disorder is displayed 
in Fig. H](b), showing well-separated quasi-modes with 
large fluctuations in p(r) along the mode. The LDOS 
is sufficient to calculate the light-matter coupling, but 
a quasi-mode description is instructive for physical in- 
sight. It is obtained by fitting the LDOS by a sum of 
Lorentzians: 



tt(uj- uj z ) 2 + (^/2) 2 



(1) 



where each term of the sum describes a single quasi 
mode with resonance frequency c^, photon decay rate 
Ki = uJi/Qi, quality factor Q^, and amplitude po,i(r). 
It is instructive to express the coupling strength in 
terms of the mode volume V{ that enters in the 
Jaynes-Cummings model [4]. For an emitter on res- 
onance with an Anderson- localized mode (uj = ujj) 
and with a transition dipole moment oriented along 
the local electric field, the Purcell factor is Fi(r) = 
67rc 3 Qpo y i(r)/uj^n(r), and from that the effective mode 
volume Vi = 1/ max r {n 2 (r)po,i(^)} is defined. The dis- 
tribution of mode volumes is subsequently calculated 
from an ensemble of 500 different disorder configurations 
for each value of the localization length and plotted in 
Fig. [2j The distribution is found to shift to smaller mode 
volumes with decreasing localization length in agreement 
with the intuitive expectation that a short localization 
length leads to small mode volume and hence efficient 
cavity QED. For a localization length of £ = 10 |Lim 
we extract mode volumes as small as Vi ~ 2.5(A/n) 3 , 
which is competitive with engineered PhC nano-cavities 
2lj. In recent experiments using quantum dots to probe 
Anderson-localized modes a localization length below 
£ = 10 |Ltm has been extracted [14| and the ultimate lower 
bound on the localization length has not yet been estab- 
lished [22|. This illustrates the very promising potential 
of efficient confinement of light with Anderson-localized 
modes. 
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FIG. 2. (color online). Histograms of the mode volume prob- 
ability distribution calculated for 500 realizations of disorder 
and for two values of the localization length £. The transverse 
extension of the modes is obtained from the simulation of the 
PhC waveguide shown in Fig. [TJa). The solid and dashed 
lines are guides to the eye. Inset: Mean values and one and 
two standard deviations of the calculated mode volume dis- 
tributions plotted as a function of localization length. 
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FIG. 3. (color online). Normalized Q-factor probability 
distributions for different localization lengths between £ = 
10 — 45 u.m and without loss (filled symbols) with fits to 
log-normal distributions (solid lines), and for a loss length 
/ = 2.5 mm and £ = 10 um (open circles) with a fit to a log- 
normal distribution modified to include loss [3 (dashed line). 
The sample length is L = 100 um. Inset: Mean value of the 
Q-factor distributions (symbols) and the central 68% of the 
Q-factor distribution (bars), with and without loss (blue dots 
and green squares, respectively). 



The Q-factor is the additional figure-of-merit deter- 
mining the light-matter coupling strength. The statistics 
of the Q-factors is obtained from an ensemble of 8000 dis- 
order configurations for each localization length leading 
to the identification of approximately 60.000 Anderson- 
localized modes near the center of the waveguide and 
within the spectral range of A = 970 — 980 nm. Fig- 
ure [3] plots the extracted Q-factor distributions for dif- 



ferent localization lengths verifying that the Q-factors of 
Anderson-localized modes in ID are described by log- 
normal distributions [H, Q (solid lines in Fig. [3]). Due 
to the long tails of the log-normal distributions, a fi- 
nite probability exists of having very large Q-factors 
that are much larger than the most probable value. 
The inset of Fig. [3] shows that the average Q-factor in- 
creases super-exponentially with an increasing localiza- 
tion length, which illustrates the very sensitive depen- 
dence of the light-matter interaction on the localization 
length that is responsible for the pronounced cavity QED 
effects observed in experiments jlQj . 

In real PhC waveguides, loss primarily due to scat- 
tering of light out of the membrane structure limits the 
quality of the light confinement. Such a loss is included 
as a finite value of the imaginary part n" of the refractive 
index n' + in" leading to the loss length / = \/2nn" and 
an associated loss Q-factor Q[ oss tt/XI, which is valid in 
the limit where loss can be treated perturbatively. This 



results in an effective Q-factor Q 



l 

eff 



' loss ' 



where 



Qioss sets the limit of the highest Q-factor that can be 
reached. An example is shown for £ = 10 /im in Fig. [3] 
(open circles) leading to a cut-off at Qioss = 28000 for 
the loss length of / = 2.5 mm. The cut-off results in an 
interesting reshaping of the Q-factor distribution since 
the highest Q-factors from the loss-less case all accumu- 
late in the part of the distribution with highest Q. The 
overall reduction in the Q-factors in the presence of loss 
is reflected in the average Q-factor that is displayed in 
the inset of Fig. [3j 

We now turn to the statistic of the light matter cou- 
pling strength and derive the probability of obtaining 
strong coupling in disordered waveguides. By apply- 
ing the decomposition into quasi- modes of Eq. (pQ), the 
threshold for reaching strong coupling for a single emitter 
that is resonant with a single Anderson- localized mode z, 
is given as 



8d 2 ' 



(2) 



The right hand side of Eq. (|2]) only contains physical con- 
stants and the dipole moment d of the emitter while the 
left hand side contains all the optical properties of the lo- 
calized modes. For identical emitters the probability for 
reaching the strong-coupling regime is therefore solely de- 
termined by the distributions of Qi and po,i that are both 
linked to the localization length, as explored in the data 
of Fig. [2] and Fig. [3l We stress that a shortening of the 
localization length is beneficial in two ways since it simul- 
taneously increases the Q-factor and decreases the mode 
volume, which is the underlying reason for the success of 
Anderson- localized cavities in cavity QED experiments. 

Based on the above criterion and the calculated dis- 
tributions we can estimate the probability that a quan- 
tum dot within a disordered PhC waveguide will be 
strongly coupled to an Anderson-localized mode. This 
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FIG. 4. (color online). Total probability of achieving strong 
coupling plotted as a function of localization length for differ- 
ent loss lengths and for a sample length of L — 100 Lim. The 
solid lines are guides to the eye. 



finite 'yield' of cavities in the strong-coupling regime orig- 
inates from the statistical nature of multiple scattering, 
but the outcome is in fact not different from the case of 
engineered cavities where fabrication imperfections need 
to be eliminated as good as possible. The resulting strong 
coupling yield in engineered cavities has to our knowl- 
edge not been studied systematically in the literature, 
but it is not unlikely that the performance of Anderson- 
localized cavities is superior to standard cavities. We 
emphasize that in the estimates below it is not assumed 
that the emitter is optimally positioned at an antinode 
of an Anderson-localized modes, which corresponds to 
the realistic experimental situation where no spatial con- 
trol of quantum dot position is implemented. We as- 
sume a transition dipole moment of d = 0.64 e nm, which 
was recently measured for InAs quantum dots [25] . The 
results are summarized in Fig. S] displaying the strong 
coupling probability as a function of localization length 
and different loss lengths for a quantum dot located in 
the center of a sample and tuned into resonance with 
the nearest Anderson- localized mode present. For a lo- 
calization length of £ > 30 |Ltm and negligible loss the 
fraction of strongly coupled systems is below 1%. This 
probability dramatically increases with decreasing local- 
ization length leading to 55% for £ = 7 |im. The rapid 
increase is mainly due to the drastic increase in the Q- 
factors while the mode volume variations are relatively 
modest. The presence of loss results in modifications 
of the Q-factor distributions (see Fig. [3j), and thereby a 
reduction of the strong-coupling probability, which is in- 
vestigated systematically in Fig. [U For £ = 7 |Ltm and 
I = 0.7 mm corresponding to the values observed experi- 
mentally in the very first experiments on disordered pho- 
tonic crystal waveguides with embedded emitters [14j], a 
strong-coupling probability of 1% is predicted. Conse- 
quently the observation of strong coupling of a single 



quantum dot to an Anderson-localized modes appears 
within experimental reach with disordered PhC waveg- 
uides presently available, and finding methods to reduce 
the localization length and increase the loss length even 
further, e.g., by inducing correlated disorder, will further 
increase this probability. In order to estimate the fun- 
damental boundaries set by out-of-plane scattering in a 
PhC waveguide we consider the simulations of Ref. [1H 
for a fixed localization length of £ = 7 |Ltm = 27a and 
realistic random perturbations in the hole diameter of 
a = 0.005a. In this case the loss length can be longer 
than I ~ 2.5 mm leading to a predicted strong coupling 
probability above 20% (Fig. \Mb)). We note that in a 
more complete model, a distribution of out-of-plane light 
leakage scattering rates is included [3, E3] as opposed to 
the single-parameter model employed here, but our con- 
clusions are expected to be robust to such extensions flBl ). 
We conclude that apart from providing a new route to 
cavity QED, Anderson localization may be competitive 
with engineered nanocavities that inherently suffer from 
the deteriorating role that disorder plays in this case. 

In summary, we have developed a theoretical model for 
a quantum dot coupled to an Anderson-localized mode. 
The probability of achieving strong coupling between a 
single photon and a quantum dot is evaluated and the 
dependence on the localization length and loss length is 
studied. Our results show that light-matter entangle- 
ment is within experimental reach in a disordered pho- 
tonic medium, and that Anderson localized random cavi- 
ties may be competitive with engineered nanocavities for 
cavity QED experiments. 
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